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Audio signal encoding or decoding 



The invention relates to encoding an audio signal or decoding an encoded 

audio signal. 

5 Erik Schuijers, Werner Oomen, Bert den Blinker and Jeroen Breebaart, 

"Advances in Parametric Coding for High-Quality Audio", Preprint 5852, 1 14th AES 
Convention, Amsterdam, The Netherlands, 22-25 March 2003 disclose a parametric coding 
scheme using an efficient parametric representation for the stereo image. Two input signals 
are merged into one mono audio signal. Perceptually relevant spatial cues are explicitly 

10 modeled as is shown in Fig. 1 . The merged signal is encoded using a mono parametric 
encoder. The stereo parameters Interchannel Intensity Difference (ED), the Interchannel 
Time Difference (ITD) and the Interchannel Cross-Correlation (ICC) are quantized, encoded 
and multiplexed into a bitstream together with the quantized and encoded mono audio signal. 
At the decoder side, the bitstream is de-multiplexed to an encoded mono signal and the stereo 

1 5 parameters. The encoded mono audio signal is decoded in order to obtain a decoded mono 
audio signal m' (see Fig. 2). From the mono time domain signal, a de-correlated signal is 
calculated using a filter D yielding perceptual de-correlation. Both the mono time domain 
signal m' and the de-correlated signal d are transformed to the frequency domain. Then the 
frequency domain stereo signal is processed with the HD, ITD and ICC parameters by 

20 scaling, phase modifications and mixing, respectively, in a parameter processing unit in order 
to obtain the decoded stereo pair P and r\ The resulting frequency domain representations 
are transformed back into the time domain. 

25 An object of the invention is to provide advantageous audio encoding or 

decoding using spatial parameters. To this end, the invention provides an encoding method, 
an audio encoder, an apparatus for transmitting or storing, a decoding method, an audio 
decoder, a reproduction apparatus and a computer program product as defined in the 
independent claims. Advantageous embodiments are defined in the dependent claims. 
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According to a first aspect of the invention, an audio signal is encoded, the 
audio signal including a first audio channel and a second audio channel, the encoding 
comprising subband filtering each of the first audio channel and the second audio channel in 
a complex modulated filterbank to provide a first plurality of subband signals for the first 
5 audio channel and a second plurality of subband signals for the second audio channel, 
downsampling each of the subband signals to provide a first plurality of downsampled 
subband signals and a second plurality of downsampled subband signals, further subband 
filtering at least one of the downsampled subband signals in a further filterbank in order to 
provide a plurality of sub-subband signals, deriving spatial parameters from the sub-subband 

10 signals and from those downsampled subband signals that are not further subband filtered, 
and deriving a single channel audio signal comprising derived subband signals derived from 
die first plurality of downsampled subband signals and the second plurality of downsampled 
subband signals. By providing a further subband filtering in a subband, the frequency 
resolution of said subband is increased. Such an increased frequency resolution has the 

1 5 advantage that it becomes possible to achieve higher audio quality (the bandwidth of a single 
sub-band signal is typically much higher than that of critical bands in the human auditory 
system) in an efficient implementation (because only a few bands have to be transformed). 
The parametric spatial coder tries to model the binaural cues, which are perceived on a non- 
uniform frequency scale, resembling the Equivalent Rectangular Bands (ERB) scale. The 

20 single channel audio signal can be derived directly from the first plurality of downsampled 
subband signals and the second plurality of downsampled subband signals. However, the 
single channel audio signal is advantageously derived from sub-subband signals for those 
downsampled subbands that are further subband filtered, in which case the sub-subband 
signals of each subband are added together to form new subband signals and wherein the 

25 single channel audio signal is derived from these new subband signals and the subbands from 
the first and second plurality of subbands that are not further filtered. 

According to another main aspect of the invention, audio decoding of an 
encoded audio signal is provided, the encoded audio signal comprising an encoded single 
channel audio signal and a set of spatial parameters, the audio decoding comprising decoding 

30 the encoded single channel audio channel to obtain a plurality of downsampled subband 

signals, further subband filtering at least one of the downsampled subband signals in a further 
filterbank in order to provide a plurality of sub-subband signals, and deriving two audio 
channels from the spatial parameters, the sub-subband signals and the downsampled subband 
signals for those subbands that are not further subband filtered. By providing a further 
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subband filtering in a subband, the frequency resolution of said subband is increased and 
consequently higher quality audio decoding can be reached. 

One of the main advantages of these aspects of the invention is that parametric 
spatial coding can be easily combined with Spectral Band Replication ("SBR") techniques. 
5 SBR is known per se from Martin Dietz, Lars Liljeryd, Kristofer Kjorling and Oliver Kunz, 
"Spectral Band Replication, a novel approach in audio coding",Preprint 5553, 112 th AES 
Convention, Munich, Germany, 10-13 May 2002, and from Per Ekstrand, "Bandwidth 
extension of audio signals by spectral band replication", Proc. 1st IEEE Benelux Workshop 
on Model based Processing and Coding of Audio (MPCA-2002), pp. 53-58, Leuven, 

10 Belgium, November 15, 2002. Further reference is made to the MPEG-4 standard ISO/EC 
14496-3 :2001/FDAM1, JTC1/SC29/WG1 1, Coding of Moving Pictures and Audio, 
Bandwidth Extension which describes an audio codec using SBR. 

SBR is based on the notion that there is typically a large correlation between 
the low and the high frequencies in an audio signal. As such, the SBR process consists of 

1 5 copying the lower part(s) of the spectrum to the higher partes) after which the spectral 

envelope is adjusted for the higher part(s) of the spectrum using little information encoded in 
the bit stream. A simplified block diagram of such an SBR enhanced decoder is shown in Fig. 
3. The bit-stream is de-multiplexed and decoded into core data (e.g. MPEG-2/4 Advanced 
Audio Coding (AAC)) and SBR data. Using the core data die signal is decoded at half the 

20 sampling frequency of the full bandwidth signal. The output of the core decoder is analyzed 
by means of a 32 bands complex (Pseudo) Quadrature Mirror Filter (QMF) bank. These 32 
bands are then extended to full bandwidth, i.e., 64 bands, in which the High Frequency (HF) 
content is generated by means of copying part(s) of the lower bands. The envelope of the 
bands for which the HF content is generated is adjusted according to the SBR data. Finally by 

25 means of a 64 bands complex QMF synthesis bank the PCM output signal is reconstructed. 

The SBR decoder as shown in Fig. 3 is a so-called dual rate decoder. This 
means that the core decoder runs at half the sampling frequency and therefore only a 32 
bands analysis QMF bank is used. Single rate decoders, where the core decoder runs at the 
full sampling frequency and the analysis QMF bank consists of 64 bands are also possible. In 

30 practice, the reconstruction is done by means of a (pseudo) complex QMF bank. Because the 
complex QMF filter bank is not critically sampled no extra provisions need to be taken in 
order to account for aliasing. Note that in the SBR decoder as disclosed by Ekstrand, the 
analysis QMF bank consists of only 32 bands, while the synthesis QMF bank consists of 64 
bands, as the core decoder runs at half die sampling frequency compared to the entire audio 
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decoder. In the corresponding encoder however, a 64 bands analysis QMF bank is used to 
cover the whole frequency range. 

Although the invention is especially advantageous for stereo audio coding, the 
invention is also of advantage to coding signals with more than two audio channels. 
5 These and other aspects of the invention are apparent from and will be 

elucidated with reference to the embodiments described hereinafter. 



In the drawings: 

10 Fig. 1 shows a block diagram of a unit for stereo parameter extraction as used 

in a Parametric Stereo ("PS 11 ) encoder; 

Fig. 2 shows a block diagram of a unit for the reconstruction of a stereo signal 
as used in a PS decoder; 

Fig. 3 shows a block diagram of a Spectral Band Replication ("SBR") decoder; 
IS Fig. 4 shows a block diagram of combined PS and SBR enhanced encoder 

according to an embodiment of the invention; 

Fig. 5 shows a block diagram of combined PS and SBR enhanced decoder 
according to an embodiment of the invention; 

Fig. 6 shows an M bands downsampled complex QMF analysis (left) and 
20 synthesis bank (right); 

Fig. 7 shows a magnitude response in dB of a prototype filter; 

Fig. 8 shows a magnitude responses in dB of the first four out of 64 non- 
downsampled complex modulated analysis filters; 

Fig. 9 shows a block diagram of a Q bands filter bank with trivial synthesis; 
25 Fig. 1 0 shows a combined magnitude response in dB of a first non- 

downsampled modulated QMF filter and 8 bands complex modulated filter bank; 

Fig. 1 1 shows a stylized magnitude response of 4 bands evenly stacked filter 
bank (top) and oddly stacked filter bank (bottom) according to an embodiment of the 
invention; 

30 Fig. 12 shows a 77 bands non-uniform hybrid analysis filter bank based on 64 

bands complex analysis QMF according to an embodiment of the invention; 

Fig. 13 shows a 71 bands non-uniform hybrid analysis filter bank based on 64 
bands complex analysis QMF for use in an audio decoder; and 
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Fig. 14 shows a block diagram of efficient implementation of the complex 
modulated analysis filter bank 

The drawings only show those elements that are necessary to understand the 

invention. 



Combining SBR with PS potentially yields an extremely powerful codec. Both 
SBR and PS are post-processing algorithms in a decoder consisting of a fairly similar 
structure, i.e., some form of time to frequency conversion, processing and finally frequency 

10 to time conversion. When combining both algorithms, it is required that both algorithms can 
run concurrently on e.g. a DSP application. Hence, it is advantageous to reuse as much as 
possible of the calculated intermediate results of one codec for the other. In the case of 
combining PS with SBR this leads to reusing the complex (Pseudo) QMF sub-band signals 
for PS processing. In a combined encoder (see Fig. 4) the stereo input signal is analyzed by 

IS means of two 64 bands analysis filter banks. Using the complex sub-band domain 

representation, a PS calculation unit estimates the stereo parameters and creates a mono (sub- 
band) down-mix is created. This mono down-mix is then fed to an SBR parameter estimation 
unit Finally the mono down-mix is converted back to the time domain by means of a 32 
bands synthesis filter bank such that it can be coded by the core decoder (core decoder needs 

20 only half the bandwidth). 

In the combined decoder as shown in Fig. 5, regardless whether or not a dual 
rate or a single rate system is being used, the full bandwidth (64 bands) subband domain 
signals after envelope adjustment are converted to a stereo set of subband domain signals 
according to the stereo parameters. These two sets of sub-band signals are finally converted 

25 to the time domain by means of the 64 bands synthesis QMF bank. If one would just combine 
PS with SBR, the bandwidth of the lower frequency bands of the QMF filter is larger than 
what is required for a high quality stereo representation. So, in order to be able to give a high 
quality representation of the stereo image, a further sub-division of the lower sub-band 
signals is performed according to advantageous embodiments of the invention. 

30 For a better understanding of aspects of the invention, the theory behind 

complex QMF sub-band filters is first explained. 
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OMF sub-band filters 

The QMF analysis sub-band filter can be described as following. Given a real 
valued linear phase prototype filter p(y) , an M -band complex modulated analysis filter 

bank can be defined by the analysis filters 

5 

h k (v) = p(v) exp + 1 /2)(y - 0) j , (1) 

for fc = 0,l> M- 1. The phase parameter 9 is not important for the analysis that follows, 
but a typical choice is (N + M )/2 , where N is the prototype filter order. Given a real 

10 valued discrete time signal x(y) , the sub-band signals v k (n) are obtained by filtering 

(convolution) x{y) with h k (y) , and then downsampling the result by a factor M (see left 

hand side of Fig. 6). 

A synthesis operation consists of first upsampling the QMF sub-band signals 
with a factor M , followed by filtering with complex modulated filters of the type (1), adding 
15 up the results and finally taking twice the real part (see right hand side of Fig. 6). Then near- 
perfect reconstruction of real valued signals can be obtained by suitable design of a real 
valued linear phase prototype filter p(y) . The magnitude response of the prototype filter as 

used in the SBR system of the MPEG-4 standard (refered to above) in case of 64 bands is 
shown in Fig. 7. Hie magnitude responses of the 64 complex modulated analysis filters are 

20 obtained by shifting the magnitude response of the prototype filter p( v) by — (k + 1 / 2) . 

M 

Part of these responses is shown in Fig. 8. Note that only the positive frequencies are filtered, 
except for k = 0 and k = M-l. As a result the sub-band signals prior to downsampling are 
close to being analytic, facilitating easy amplitude and phase modifications of real-valued 
sinusoids. Phase modifications are also possible for the first and last band as long as the 
25 sinusoids residing in these bands have a frequency that is above n 1 2Af or below n-nl 2M 
respectively. For frequencies outside this region the performance of phase modification 
deteriorates rapidly because of interference of the negative frequencies. 

Starting from the QMF analysis filters as described above, in embodiments of 
30 the invention, a finer frequency resolution is obtained by further filtering each downsampied 
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subband signal v k (w) into Q k sub-subbands. In the following the properties of the further 
subband filtering will be derived. 

Signal modification in the complex OMF sub-band domain 
5 In the following, let Z(a>) = 2] z ( n ) exp(-/>wy) be the discrete time Fourier 

transform of a discrete time signal . Assuming the near-perfect reconstruction property 
as mentioned above and also a design where P{a>) , the Fourier transform of p(y) , 
essentially vanishes outside the frequency interval \-nlM>nlM\ , which is the case for the 
prototype filter p(y) as illustrated above, the next step here is to consider a system where the 
1 0 sub-band signals v k (n) are modified prior to synthesis. Now, let each sub band k be 
modified by filtering with a filter B k (aj) . With the extending definition 

B*(fi>) = for *<0, (2) 

where the star denotes complex conjugation, it can then be shown (neglecting overall delay, 
assuming a real valued input and a single rate system) that the resulting system including 
1 5 filter bank synthesis corresponds to a filtering with the filter 

B{&) = B k {Ma>) \P(c> - n(k + 1/2)/ M)f . (3) 

According to the hypotheses regarding the properties of P(o?) , inserting 
B k {<o) = 1 for all k in (3) leads to B(<d) = 1 , and a squared sum identity follows for the 
shifted prototype filter responses. By choosing real-valued constants B k (jCO) = b k £ 0 the 
20 system acts as an equalizer, which interpolates the gain values b k at frequencies 

n{k + 1 / 2) / M . The attractive feature is that the overall system is time-invariant, that is, free 

of aliasing, in spite of the use of down- and upsampling. This will of course only be true up 
to the amount of deviation to the stated prototype filter hypotheses. 

In order to derive a mono audio signal, additional sub-filtering of the complex 
25 sub-band signals should not only preserve these properties, but also extend these properties to 
manipulation of the filtered sub-band signals. Sub-filtering preserving these properties can be 
performed using a modification of so-called Mth band filters as known per se from P.P. 
Vaidyanathan, "Multirate systems and filter banks", Prentice Hall Signal Processing Series, 
1993, sections 4.6.1-4.6.2). 
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Modulated filter banks with trivial synthesis 

A discrete time signal v(n) can be split into Q different signals by a bank of 

filters with impulse responses g q (n) 9 q = 0,l,...,fi-l - This is illustrated in Fig. 9. 
Let the corresponding analysis outputs be y q (ji) , and consider the trivial synthesis operation 

yOO = • (4) 

Perfect reconstruction, y(ri) = v(n) , is then obtained by choosing the filters such that 

4 00 -'00» (5) 

where S(n) = 1 if n = 0, and S(n) = 0 if n * 0 . For causal filters, the right hand side of (5) 
would have to be replaced with S (n - d) where d is a positive delay, but this straightforward 
1 0 modification is omitted for clarity of exposition. 

The filters g (/?) can be chosen as complex modulations of a prototype filter g(n) through 




g q (n) = *(n)expjf^(g + l/2)nj . (6) 

In this preferred embodiment of the invention, the filters are oddly stacked (the 
1 5 factor q + l/2). An advantage of this preferred embodiment will be explained later. Perfect 

reconstruction (S) is obtained if and only if 

g(Qn) = S(n)/Q. (7) 
A variation of this is the real-valued cosine modulation as 

g q <ji) = £(*)cosj^(9 + l/2)*j, (8) 

20 with a real-valued prototype filter g(m) satisfying 

g(2Qn) = S(n)/Q. (9) 
(This is easily obtained by consideration of g (n) + g Q ^ (n) in (6).) 



Sub-filtering the complex-exponential modulated filter bank 
25 Starting from the QMF analysis filters as described above, a finer frequency 

resolution is obtained by further filtering each downsampled subband signal v k (/)) into Q k 

sub-subbands by using one of the modulated structures (6) or (8) above. Denote the resulting 
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output signals >/* (/*) , and let g k (n) describe the filter bank applied within sub band k . If 

Q k = 1 , there is no filtering and g^i/i) = 5(ri) . A typical application example is the case 
where M = 64, & =8, Q k = 4 for * = 1,2, and Q k =lfor k>2. 

The combined effect of the two filter banks from x(y) to y k (n) can be 
5 described as filtering with filters F k ((D) followed by downsampling by a factor M , where 

F k {a>) = H k (a>)G k q (Ma>) . (10) 
If the prototype filter response P(&) is essentially zero outside the interval \-x I M y nl M ] , 
which is the case for the SBR analysis filters (see Fig. 7), then the filter F k (o) has a single 

nominal center frequency defined in the complex modulated case by 
10 ofc, =27v(q + Q k s+$/(MQ k ) 9 (11) 

where s is a integer chosen such that Q k {k — £) £ 2(? + g^) + 1 £ +•£) . For example, as 
illustrated in Fig. 10, if k = 0 and Q Q = 8 , the values of q? 00 ,(D 01> ... 9 q> 01 are 

TV 



(1,3,5,7,9,11,-3,-1). 



15 Signal modification with non-uniform frequency resolution 

The insertion of sub-subband filter banks as described above does not 
introduce further downsampling, so the alias-free performance of signal modification as 
shown above in the case of complex QMF only, is preserved. Consider the general combined 
operation of M -subband analysis, further subband filtering by using Q k sub-subbands 

20 within subband k , filtering of each sub-subband signal y k (n) by a filter ^ (a>) , synthesis 

within each subband k by summation, and finally synthesis through the M -band synthesis 
bank. The overall transfer function of such a system is given by (3) with, for k > 0 , 

B k {*>) = X AA*>)G k {co) , (12) 

For a> > #7(2 A/) , this gives 

25 B(^)= £^^(M*>) G k (Ma))\P(a>-7r{k + l/2)/M)\ 2 9 (13) 

so the throughput response of the sub-subband (*, q) is G q {Mm) \P(& - n{Jk + 1 / 2) / M)| 2 . 
For \a\ < rrf(2M) , some care has to be taken due to (2). In this frequency range it holds that 
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B(a>) = B Q (Ma>)\P(<D - * /(2M))f + B q (-Mq>Y |P( /(2M))f (14) 



and assuming a real sub-subband prototype filter coefficients, it holds that 



G° q (-a>y =Gl+ q (a>) % (15) 



so if the modifying filters are chosen such that 




then B 0 (-Ma)Y = B 0 (Mw) and the squared sum identity mentioned in connection with (3) 
leads to 



and negative frequencies. This is the reason why oddly stacked (complex) filters are being 
used for sub-filtering the QMF sub-band signals instead of evenly stacked (complex) filters 
(see Fig. 11). For evenly stacked filters it is not possible to apply phase modifications of 
sinusoids residing in the centre filter, i.e., the filter with a centre frequency of zero, as there is 
15 no discrimination between positive and negative frequencies possible. Assuming a prototype 
filter with a response G{a>) band limited to [-2fl7g,2;r/g], with Q the number of bands, 

for the evenly stacked case the lower limit to which phase modifications can approximately 
be applied is Ix/Q, whereas for the oddly stacked case the lower limit to which phase 

modifications approximately can be applied is n I Q . 

20 As mentioned in the introduction, for PS synthesis important special cases of 

the above are equalization and phase modification For equalization, 
4*0°) = a k4 ^ 0 and the condition (16) reduces to 




(17) 



10 



for|o)| 7r/(2M) , corresponding to a throughput response G*(Ma>) for sub-subband (0,^) . 

Equations (15) until (17) indicate the desire to discriminate between positive 



(18) 




30 



Stereo parameter estimation 

The non-uniform complex filter bank, i.e. the QMF bank followed by the 
further subband filtering, as described above, can be applied to estimate the stereo parameters 
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Inter-channel Intensity Differences (OD), Inter-channel Phase Differences (DPD) and Inter- 
channel Cross Correlation (ICC) as shown below. Note that in this practical embodiment, 
IPD is used as a practically equivalent substitute for the 11 D as used in the paper of Schuijers 
et al. In the combined PS encoder (see Fig. 4) the first three complex QMF channels are sub- 
5 filtered so that in total 77 complex-valued signals are obtained (see Fig. 12). 

From this point on the 77 complex-valued time-aligned left and right sub- 

subband signals are denoted as /*(*) and r*(«) respectively, accordingly the indexing of 

To estimate the stereo parameters at a certain sub-band sample position ri the 
10 left, right and non-normalized cross-channel excitation are calculated as: 




for every stereo bin b , h(n) is the sub-band domain window with length L , e. a very small 
value preventing division by zero (e.g. s = le - 10 ) and /* (n) and r* (n) the left and right 
sub-subband domain signals. In case of 20 stereo bins, the summation over k from k t up to 
1 5 and including k h and q from q { up to and including q h goes as shown in Table. Note that 
the 'negative* frequencies (e.g. £ = 0 with q - 4.. .7) are not included in the parameter 
estimation of (20). 
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Table 1: Start and stop indices of summation over k and q 
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19 
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63 


0 


0 


35*/64-* 



The summations to calculate £/(&), and are aligned such that the mid-point of 
these signals in the summation coincides with the parameter position, hence the shift by 

As is clear from Table 1, only sub-subband signals and subband signals with a 
positive centre frequency are used for estimating stereo parameters. 



The IE), denoted as l(b), the ICC, denoted as C(b) and the IPD, denoted as P(b) for each 
stereo bin b are calculated as: 
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'(»)=101o gl0 (f$>) 



(21) 



P(b)=^e R {b) 

The angle in the equation P(b) = Ze R (b) is calculated using the four quadrant arctangent 

function giving values between - n and n . Depending on target bit rate and application, 
these parameters, or a subset of these parameters are quantized and coded into the PS part of 
5 the bit-stream. 



10 



15 



20 



Stereo signal synthesis 

In order to keep the computational costs (in terms of RAM usage) in the 
decoder as low as possible a similar analysis structure is used. However the first band is only 
partially complex (see Fig. 1 3). This is obtained by summation of the middle band pairs 

G£ (<y) and G° (<*>) and G 3 °(a>) and G° {col) . Furthermore, the second and the third band are 
two-band real-valued filter banks, which is obtained by summation of the output of G * (a>) 

and G* (co) , and summation of the output of G* (&) and G* (a?) (see also the discussion in 

the section about modulated filter banks). Using this simplification of the decoder filter-bank 
structure, still the discriminative feature between positive and negative frequencies is 
maintained by subdivision of the first sub-band filter. The decoder analysis filter bank is 
shown in Fig. 13. Notice that the indexing of the first QMF filtered (sub-)subband signals is . 
sorted according to frequency. 

The stereo (sub-)subband signals of a single frame are constructed as: 

^W=A n ^(n)+A 21 £/ A (/i) 
r k (n) = A 12 ^ (n)+ (n) 



(22) 



(23) 



with s k (n) the mono (sub-)subband signals, and d k (n) the mono de-correlated (sub- 
)subband signals that are derived from the mono (sub-)subband signals s k (ri) in order to 
account for synthesizing the ICC parameters, k = 0,..., K - 1 the sub-band index (K is the 
25 total number of sub-bands, i.e., K = 71), QMF sub-band sample index n = 0,..., Af-1 with N 
the number of sub-band samples of a frame, A n , A l2 , A 21 , A 2 the scale factor 



WO 2005/043511 PCT/IB2004/052226 

14 



5 



manipulation matrices and P„ the phase rotation manipulation matrix. The manipulation 

matrices are defined as function of time and frequency and can be derived straightforwardly 
from the manipulation vectors as described in the MPEG-4 standard ISO/IEC 14496- 
3:2001/FPDAM2, JTC1/SC29/WG1 1, Coding of Moving Pictures and Audio, Extension 2. 

s k (n) is defined according to Fig. 12 as resulting in Fig. 13: 

*o 00 ->•(») 
*i(»)-J?0») 
»a0O->JO0 
*s 00 

^(»)=v J 0 («)+v 5 °(n) 

'sOO-jJOO+^OO (24) 

*6 00 = ^ 00 + ^1 00 

*» 00 -*}(*)+>£ 00 

*•(») = >2 00 +>? 00 
s 9 (n) = yf(n)+yl(n) 

'tOO-jT'OO k = 10.-70 

Synthesis of the stereo parameters takes place accordingly the indexing of Table 1. 



10 Table 1: Parameter indexing table 
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9 


7 


5*/128-6*/128 


10 


8 


3*/64-4*/64 


11 


9 


4^/64-5^-/64 
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19 


10 

lv 


/ AA — 6*r / AA 
DTI / \>*t U/fr / Q*f 


ID 


1 1 

1 1 


/AA — 7fr / AA 
O/l 1 0*t — 111 1 OH 


1 A 


19 


7*r / AA — 8*r / AA \ 


1 D 






1 £-17 


1 A 


Qfr / AA _ 1 W / AA 
/ OH — 1 1/t / OH 




ID 


1 1 rr / 64 — 1 4tt / 64 


21-24 


16 


14^/64-18^/64 


25-29 


17 


18*764-23/r/64 


30-41 


18 


23*7 64 -35*/ 64 


42-70 


19 


35«764-;r 



(25) 



The synthesis equations thus look like: 

'*(») = A u OW,nK(n)+A 2 ,(/W,nK(«) 
'*(»)= A 12 (/(*),wK(»)+ A a (i(*),iiK 00 

/*(«) = '*(«> yP - ( ' ( * V,) 
'•*(") = '*(»>" yP ** ( ' ( * X " ) 

Note that the sign of P rt changes in the equations above if a * is encountered in the table. 

This is accordingly equation (19), i.e., the inverse phase rotation has to be applied for the 
negative frequencies. 



(26) 



Efficient implementation of modulated filter banks with trivial synthesis 
1 0 Given a modulated filter bank with a prototype filter of length L , a direct 

form implementation would require QL operations per input sample, but the fact that the 

modulation in (6) is antiperiodic with period Q can be used to split the filtering into a 

polyphase windowing of L operations followed by a transform of size Q for each input 

sample. Please note that a polyphase representation as such is known from P.P. 
15 Vaidyanathan, "Multirate systems and filter banks", Prentice Hall Signal Processing Series, 
1993, section 4.3). The following provides an advantagous application of such a polyphase 
representation according to a preferred embodiment of the invention. 

The transform is a DFT followed by a phase twiddle, which is of the order of 
Q 1°S2 Q > w hen Q is a power of two. So a large saving is obtained in typical cases where 
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L is much larger than log 2 Q . In the real modulated case (8), antiperiodicity of period 

2Q combined with even/odd symmetries around « = 0 and n = Q can again be used for 

polyphase windowing, and the transform kernel is a DCT of type HI. A detailed description 
for the case of complex modulation is given below. 
5 An effective implementation of the sub-subfiltering, using FFT core 

processing, may be realized using poly-phase decomposition of the prototype filter followed 
by modulation. Assume a prototype filter g(ri) of order N , where N = mQ and m is a 

positive integer. This condition is not restrictive, since a prototype filter of arbitrary order can 
be zero padded to fulfill the constraint. The Z-transfonn of the prototype filter designed for 
10 use in a complex modulated system (6) is 



N/2 



G{z)= 2 (27) 



This may be expressed in poly-phase notation as 



GOO = (28) 



where 



15 £,(z)= *ff g(Qn+l)z- n (29) 



All filters of the filterbank are frequency-modulated versions of the prototype filter. The Z- 
transformof the filter g q {ri) is given by 



GAz) = G(zW q ^) (30) 



where 



20 W = e Q (31) 

The expression for the output from one filter is 

W = W00 = G(zW^)V(z) = 

= ^E l ( r z Q )V{z)z-'W~ ri W-''' = (32) 

lot) 

tl j»t i2*»/ 

]E l (-z Q )V{z)z-'e^ e a 

By identifying the components of the last sum, it may be seen that the poly-phase 
25 components process delayed versions of the input signal, which subsequently are multiplied 
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by a complex exponential. Finally, all the output signals Y q (z\ q = 0 0-1 , are found by 
applying an inverse FFT (without scaling factor). Fig. 14 shows the layout for the analysis 
filter bank. Since the poly-phase filters in (29) are non-causal, a proper amount of delay has 
to be added to all the poly-phase components. 
S It should be noted that the above-mentioned embodiments illustrate rather than 

limit the invention, and that those skilled in the art will be able to design many alternative 
embodiments without departing from the scope of the appended claims. In the claims, any 
reference signs placed between parentheses shall not be construed as limiting the claim. The 
word 'comprising 9 does not exclude the presence of other elements or steps than those listed 

10 in a claim. The invention can be implemented by means of hardware comprising several 
distinct elements, and by means of a suitably programmed computer. In a device claim 
enumerating several means, several of these means can be embodied by one and the same 
item of hardware. The mere fact that certain measures are recited in mutually different 
dependent claims does not indicate that a combination of these measures cannot be used to 

15 advantage. 



